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1. Introduction 



The QCD phase diagram is a matter of intense investigation from both, theory and experiment. 
Some of the pressing open questions discussed at this conference are the presence or absence of 
a critical point the possibility of a confined chirally symmetric ('quarkyonic') phase and 
the (non-)coincidence of the chiral and the deconfinement transition at zero chemical potential 
Q]. Answers to these questions certainly require nonperturbative approaches to QCD. On the 
other hand, the considerable complexity of these questions suggests that one approach alone is 
hardly capable to provide all answers. Instead it seems promising to combine the various available 
methods in order to balance their respective strengths and weaknesses. 

Lattice Monte-Carlo simulations are well behaved at zero or imaginary chemical potential, 
but encounter the notorious sign problem when it comes to real chemical potential. One approach 
to overcome this problem, extrapolation from zero chemical potential by Taylor-expansion meth- 
ods (see e.g. [Q]), has been questioned recently on the basis of the failure of corresponding ex- 
trapolations in model calculations [|6|]. These models, notably the Polyakov-NJL model and the 
Polyakov-quark-meson model (see e.g. [^, ^, ^ and refs. therein) have been employed frequently 
to explore the details of the QCD phase diagram at zero and finite chemical potential. Their success 
is demonstrated e.g. by the quantitative reproduction of the lattice equation of state at zero chem- 
ical potential. Furthermore they serve as a formidable qualitative playground to explore scenarios 
for the details of the interplay between the chiral and the deconfinement transition. Nevertheless, it 
is hard to see how the various model parameters can be constrained enough to arrive at quantitative 



predictions as e.g. the location of a possible critical point [10]. 



A third class of approaches are functional methods involving the renormalization group equa- 
tions [11 1 and/or Dyson-Schwinger equations [ p^ , |l3| ] of QCD. In the past years much progress 
has been made in extending these methods to finite temperature and/or chemical potential, see e.g. 



[ ]14[ , |12| , |15| , |16| , \n[ |18| , |19| , gOp. Most of the earlier works involving functional methods concen- 
trated on the chiral aspects of the QCD transition. Recently, methods became available that also 
take into account the deconfining aspect [14, 16, 18|, 1^, 20 1. In this talk we report on results from 
one particular method that extracts the deconfinement transition temperature from the properties of 
the quark propagator at generalized boundary conditions. This method has been introduced origi- 



nally within the lattice framework [21, 22] and adapted to functional methods in [18, 19, 20]. The 
quantity signaling the deconfinement transition is the dual quark condensate (or 'dressed Polyakov 
loop'). It transforms under center transformations exactly like the ordinary Polyakov loop and is 
therefore an order parameter in the limit of infinitely heavy quarks. This quantity is furthermore in- 
teresting from a formal perspective since it relates both the chiral and the deconfinement transition 



to the spectral properties of the Dirac operator [^ 



In the following we first recall the defining equations for the ordinary and the dual quark 
condensate, then summarize the truncation scheme used in our DSE calculations before we discuss 
our results for the chiral and deconfinement phase transition. 
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2. The dual quark condensate 

Consider the quark propagator S{p, (Op) at finite temperature given by tlie tensor decomposition 

S{p,(Op) = [iY4(OpC{p,(Op)+iYiPiA{p,C0p) +B{p,(Op)r^ , (2.1) 

with vector and scalar quark dressing functions C,A,B. At physical, antiperiodic boundary con- 
ditions the corresponding Matsubara frequencies are given by (Op{nt) = {2nT){nt + 1/2). The 
ordinary quark condensate can be extr acted from the trace of the quark propagator by 

(W)^ =Z,„Z2A^cr£ j -^^tVDS{p,(Op) . (2.2) 

For vanishing bare quark masses this integral is well-behaved and delivers the chiral condensate, 
whereas at finite bare quark masses it is quadratically divergent and needs to be properly regular- 
ized. 

Consider now non-standard, ?7(l)-valued boundary conditions in the temporal direction es- 
tablished by the equation 1/r) = e"''vA(x,0) for the quark field i/a with the boundary angle 
9 G [0,27r[. For the physical antiperiodic fermion boundary conditions we have cp = n, whereas 
9 = corresponds to periodic boundary conditions. The corresponding Matsubara frequencies are 
given by (Op{n,,(p) = {2nT){nt + (p/2n). In this non-standard framework one can evaluate a quark 
condensate by 

{W)<p=Z2NcTY, J -^tr DS{p,(Op{(p)) (2.3) 



n, 



with the conventional quark condensate obtained for 9 = tt and multiplication with Z,„. It has been 



shown in [22] that the Fourier-transform 



^i=£'^e-''P{Wv)<p (2.4) 

of this 9-dependent condensate delivers a quantity that transforms under center transformations 
exactly like the Polyakov-loop and is therefore an order parameter for the deconfinement transition. 
This quantity Zi is called the dual condensate or dressed Polyakov loop. 

The relation of £1 to the ordinary Polyakov-loop becomes transparent in the following loop 
expansion of the 9-dependent condensate: 

i(pn{l) 

Here ^ denotes the set of all closed loops I with length |/| on a lattice with lattice spacing a. Fur- 
thermore m is the quark mass. U (/) stands for the chain of link variables in a loop / multiplied with 



appropriate sign and normalization factors, see Ref. [ [22| ] for details. Each loop that closes around 



the temporal boundary picks up factors of e^"'' according to its winding number «(/). Correspond- 



ingly, the Fourier-transform in Eq. ( |2.4| ) projects out exactly those loops which wind once around 
the temporal direction of the lattice (therefore the notation Zi). In the limit of heavy quark masses 
long loops are suppressed by 1 /mi'l and therefore only the straight line along the temporal direction 
of the lattice survives; the dual condensate is then equal to the ordinary Polyakov-loop. 

An interesting property of the dual condensate Zi is the fact that it can be evaluated with 
functional methods [|T^], as we will see in the following. 
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Figure 1: The Dyson-Schwinger equation for the quark propagator. Filled circles denote dressed propaga- 
tors whereas the empty circle stands for the dressed quark-gluon vertex. 



3. The Dyson-Schwinger equation for the quark propagator at finite temperature 

The Dyson-Schwinger equation for the quark propagator Eq. ( ^ ) is displayed diagrammati- 
cally in Fig. [I[ At finite temperature T it is given by 

S-\p)=Z2S,\p)-CFZifg^TY^ J ^Y^S{k)ry{k,p)D^,{p-k), (3.1) 



with p = {p,COp) and k = (k,(Ok) and renormalization factors and Zif. Here D^y denotes the 
(transverse) gluon propagator in Landau gauge and Fy the quark-gluon vertex. The bare quark 
propagator is given by Sq^{p) = iy-p + m. The Casimir factor Cp = {N^ — ^)/Nc stems from the 
color trace; here we only consider the gauge group SU{2). The quark dressing functions A,B,C 
can be extracted from Eq. ( p| ) by suitable projections in Dirac-space. 

In order to solve this equation we have to specify explicit expressions for the gluon propagator 
and the quark-gluon vertex. At finite temperatures the tensor structure of the gluon propagator 
contains two parts, one transversal and one longitudinal to the heat bath. The propagator is then 
given by {q = {q,(Oq)) 

D,.{q) = ^PlM + ^PllM (3.2) 
q H 

with transverse and longitudinal projectors 

Plviq) = (5,7 - ^) 5i^5jy , Pli,{q) = P^y{q) - PlM) , (3-3) 

with = 1 ... 3). The transverse dressing Zriq, COq) is also known as magnetic dressing function 
of the gluon, whereas the longitudinal component Ziiq) is called electric dressing function of the 
gluon propagator. At zero temperatures Euclidean C?(4)-invariance requires both dressing functions 
to agree, i.e. Zriq) = Ziiq) = Z{q). 

The temperature dependence of the gluon propagator can be inferred from recent lattice calcu- 
lations. The results of Ref. 1 24] are shown in Fig. ^ The temperature effects on both the magnetic 



and electric dressing functions are such that there are almost no effects when comparing the T = 
result with T = 119 MeV. Further increasing the temperature to T = 298 MeV and T = 597 MeV 
significantly decreases the bump in the magnetic dressing function around p^ = I GeV^. There is 
no indication that this decrease takes special notice of the critical temperature ~ 300 MeV for 
quenched QCD with gauge group SU (2). The opposite seems to be true for the electric part of the 
propagator. Here from T = 119 MeV to T = 300 MeV one observes a clear increase of the bump 
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transverse gluon dressing function 
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Figure 2: Quenched SU (2) lattice resuhs [g^ for the transverse dressing function Zj- (q) and the longitudinal 
dressing function Zi{q) of the gluon propagator together with the fit functions [jH 



in the dressing function Zi(^) and a subsequent decrease when the temperature is further raised to 
T = 597 MeV. Pending further investigation it seems reasonable to assume that the maximum of 
the bump is reached at or around the critical temperature ~ 300 MeV. 



Although the lattice data still have considerable systematic errors [ ]24| ] they may very well cor- 
rectly represent the qualitative temperature dependence of the gluon propagator. We therefore use 
a temperature dependent (qualitative) fit to the data as input into the DSE; this fit is also displayed 
in Fig. ^ (straight lines). The fit functions are described in detail in Refs. [18, 19] and shall not be 
repeated here for brevity. Note, however, that we also inherit the scale determined on the lattice 
using the string tension ^/a = 0.44 GeV [24]. 

For the quark-gluon vertex with gluon momentum q = {q,(Oq) and the quark momenta p = 
(p, (Op),k = (k, (Ok) we employ the following temperature dependent model 

, , X ^ (s, C{k)+C[p) , „ A{k)+A{p) 



d\ 



+ ■ 



q^ /A)«(At)ln[^VA2 + l] 



d2 + q^ A^+q 



471 



25 N 



(3.4) 



where 5 = —9/44 is the anomalous dimension of the vertex. The dependence of the vertex on 
the quark dressing functions A and C is motivated by the Slavnov-Taylor identity for the vertex. 
The remaining fit function is purely phenomenological, see e.g. [ p5| ] where an elaborate version 
of such an ansatz has been used to describe meson observables. The parameters are given by 
di = 7.6GeV^ and d2 = 0.5 GeV^. A moderate variation of these parameters shifts the critical 
temperatures of both, the chiral and the deconfinement transition but leaves all qualitative aspects 
of the results presented below unchanged. 

The truncation scheme described above has the merit to explicitly implement a realistic tem- 
perature dependence of the gluon propagator and the quark-gluon vertex beyond simple ansaetze, 
see e.g. [14, |12[, ^] for previous approaches. The explicit expressions of the resulting DSEs 
for the quark dressing functions together with the details of our numerical method are given in 
Ref. [j^j. 
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T ^ 1 ^ , ^ , ^ , ^ , ^ , ^ r 




T [MeV) T [MeV] 

Figure 3: Left diagram: Temperature dependence of the dressed Polyakov-loop Ei and the conventional 
quark condensate Aj^ = (xj/y/) (p=ji together with their derivatives for m ~ lOMeV. Right diagram: The same 
quantities in the chiral Hmit. 



4. Numerical results 

In Fig. |3| we display our numerical results for the ordinary and the dual quark condensate 
together with their (normalized) temperature derivatives once evaluated for a quark mass of m = 
lOMeV and once evaluated in the chiral limit. One clearly sees the difference in the chiral tran- 
sition: whereas at finite bare quark mass we encounter a crossover the transition changes into a 
second order phase transition in the chiral limit. In the first case the corresponding temperature 
derivative shows apeak at = 301(2) MeV, whereas it diverges at Tc = 298(1) MeV in the second 
case. We also extracted the corresponding transition temperatures from the chiral susceptibility 

XR = m^^({YY)T-{W)T=o) ■ (4.1) 
dm \ J 

The results for quark mass m = 10 MeV are given in table |l|. 

The corresponding transition temperature for the deconfinement transition can be read off the 
dual quark condensate (or dressed Polyakov loop). At finite quark mass and in the chiral limit we 
observe a distinct rise in the dual condensate around T w 300 MeV. The corresponding (normalized) 
temperature derivative shows peaks at T^^c = 308(2) MeV for quark mass m = 10 MeV. In the chiral 
limit this peak moves to T^ec = 299(3) MeV. 

In general we note that the chiral and deconfinement transition are close together. There are a 
few MeV difference between the different transition temperatures for the crossover at finite quark 
masses, whereas both transitions occur at the same temperature (within errors) in the chiral limit. 



These findings agree with early expectations from lattice simulations [28]. 



Tc 






1 Tiigc 


301(2) 


304(1) 


305(1) 


308(2) 



Table 1: Transition temperatures for the chiral and deconfinement transition for quark mass m — lOMeV. 
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Figure 4: Angular dependence of the quark condensate evaluated at two different quark masses and in the 
chiral hmit at T = 400 MeV. 



Furthermore we wish to emphasize that the present calculation, although carried out with 
quenched lattice results for the gluon propagator, is in itself not strictly quenched: our ansatz 
for the quark-gluon vertex is too simple to strictly represent the quenched theory. This can be seen 
from the fact that the dressed Polyakov-loop is not strictly zero below the deconfinement transition. 
Consequently we do not observe the second order deconfinement phase transition expected from 
quenched SU (2) Yang-Mills theory. Note, however, that even if our vertex were strictly quenched 
it is not clear whether the lattice input for the gluon propagator is precise enough to allow for an 
observation of the second order phase transition. 

The details of the mass dependence of the (p-dependent condensate are studied in Fig. Q We 
compare the angular dependence of the condensate at T = 400 MeV for two different quark masses 
and in the chiral limit. We clearly see a broadening in the central dip of the graphs with decreasing 
quark mass. This can be readily understood from the loop expansion of the quark condensate. 



Eq.(2.5). At sufficiently large quark masses large loops are suppressed by powers of 1/m. As a 



result only loops winding once around the torus should contribute in Eq. (2.5) and the resulting 
angular behavior of the condensate should be proportional to cos(<p). Indeed, this is what we see: 
the result for our largest quark mass can be well fitted by only few terms in an expansion A{(p) = 
lln=o^ncos{n(p) and the first term is by far the largest contribution. For smaller quark masses we 
observe also sizeable contributions from terms cos(«(p) with « > 1. In the plot, these contributions 
are responsible for the flat area around the antiperiodic boundary angle cp = n. Approaching the 
chiral limit this area becomes flatter and finally develops a derivative discontinuity at two finite 
values of (p = n ± L. These indicate the breakdown of the loop expansion Eq. ( ^!^ ) in the chiral 
limit. 

Finally we show the angular and temperature dependence of the (p-dependent condensate 
A^(r) in Fig. |[ The 3d-plot clearly shows the different evolution of the condensate at varying 
boundary conditions. Whereas at physical antiperiodic boundary angle <p = tt we observe the sec- 
ond order chiral phase transition also shown in Fig. ^, we find a monotonically rising condensate 
at the periodic boundary conditions (p = 0. For larger temperatures (not shown in the plot) we can 
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Figure 5: A 3d-plot of the angular and temperature dependence of the chiral quark condensate. 



extract a quadratic rise of the 9-dependent condensate A(p=o{T), 

A^=q{T) ~ for r > r. 



(4.2) 



This behavior can also be extracted analytically from Eqs. ( |3.1| ) and ( |2.2| ) for the quark propagator 
and the quark condensate as shown in the appendix of Ref. Around the physical value of 

(p = 71 we see a plateau with A(p{T) = that gets broader with increasing temperature. The width 
of this plateau seems to settle at a finite value smaller than In for T > 2Tc; however from the 
available results we can neither show nor exclude that it approaches In very slowly for T 00. 



5. Summary 

In this talk we addressed the chiral and the deconfinement transition of quenched QCD. We 
showed results for the order parameter for the chiral transition, the quark condensate, and an order 
parameter for the deconfinement transition, the dressed Polyakov loop extracted from the Landau 
gauge quark propagator evaluated at a continuous range of boundary conditions for the quark fields. 
We found almost coinciding transition temperatures for the chiral and the deconfinement transition 
at a moderate quark mass of the order of an up-quark. In the chiral limit the two transitions co- 
incide within error. We find a second order chiral phase transition at T^^jja = 298(1) MeV and a 
similar temperature for the deconfinement transition, T^^^ = 299(3) MeV. It is worth to emphasize 
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again that both transition temperatures are extracted from the properties of the quark propagator, 
respectively the underlying properties of the Dirac operator. 

The framework used in this work is quenched SU (2) Yang-Mills theory. Our transition tem- 
perature may be translated into the corresponding ones of quenched SU (3) QCD using the relations 
Tc/ \/o = 0.709 {SU (2)) and T^/ \/g = 0.646 {SU (3)) between the respective critical temperatures 
and the string tension [0O|]. The resulting transition temperature is then T^^jj^ « T^ec ~ 272 MeV 
in the chiral limit. In order to work in the full, unquenched theory we would have to take into ac- 
count quark-loop effects in the gluon propagator and meson effects in the quark-gluon vertex [^]. 
These effects will shift the transition temperatures below T = 200 MeV, see [||, for latest results 
fox Nf = 2+\ quark flavors. As concerns the dual condensate and scalar dressing function in the 
unquenched formulation one needs to carefully take into account effects due to the Roberge-Weiss 
symmetry p9|]. This is because of the formal similarity of the continuous boundary conditions for 
the quark field to an imaginary chemical potential, see [ 2C ] for details. 

Note: After completion of this proceedings contribution also first results for the dual conden- 
sate evaluated in the PNJL-model became available [31]. 
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